Universal Heat Conduction in the Iron- Arsenide Superconductor KFe 2 As 2 

Evidence of a rf-wave State 
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The thermal conductivity k of the iron-arsenide superconductor KFe2As2 was measured down 
to 50 mK for a heat current parallel and perpendicular to the tetragonal c axis. A residual linear 
term at T — > 0, Ko/T, is observed for both current directions, confirming the presence of nodes 
in the superconducting gap. Our value of kq/T in the plane is equal to that reported by Dong et 
al. [Phys. Rev. Lett. 104, 087005 (2010)] for a sample whose residual resistivity po was ten times 
larger. This independence of ko/T on impurity scattering is the signature of universal heat transport, 
a property of superconducting states with symmetry-imposed line nodes. This argues against an 
s-wave state with accidental nodes. It favors instead a d-wave state, an assignment consistent with 
five additional properties: the magnitude of the critical scattering rate F c for suppressing T c to 
zero; the magnitude of Ko/T, and its dependence on current direction and on magnetic field; the 
temperature dependence of k(T). 



PACS numbers: 74.25.Fy, 74.20. Rp,74.70.Dd 

The pairing mechanism in a superconductor is inti- 
mately related to the pairing symmetry, which in turn 
is related to the gap structure A(k). In a e?-wave state 
with d x 2__ y 2 symmetry, the order parameter changes sign 
with angle in the x-y plane, forcing the gap to go to zero 
along diagonal directions (±k x —±k y ). Those zeros (or 
nodes) in the gap are imposed by symmetry. The gap 
in states with s-wave symmetry will in general not have 
nodes, although accidental nodes can occur depending on 
the anisotropy of the pairing interaction. In iron-based 
superconductors, the gap shows nodes in some materials, 
as in BaFe 2 (Asi_ a; P 2 ;)2 [I] and Ba(Fei_ a; Ru a ;)2As2 
and not in others, as in Bai_ x K a; Fe2As2 [HI H] and 
Ba(Fei_ x Coa;)2As2 [5J [6] at optimal doping. 

In KFe 2 As 2 , the end-member of the Bai_ ;E K. I ;Fe 2 As2 
series (with x = 1), the presence of nodes was detected by 
thermal conductivity [7J , penetration depth [5] and NMR 
[10] . The question is whether those nodes are imposed 
by symmetry or accidental. Calculations differ in their 
predictions [TTltTB"] . Some favor a d-wave state [H], oth- 
ers an s-wave state with accidental line nodes that run 
either parallel to the c axis [TS] or perpendicular [TTJ. 
One can distinguish a cZ-wave state from an extended s- 
wave state with accidental nodes by looking at the effect 
of impurity scattering [JB]. Nodes are robust in the for- 
mer, but impurity scattering will eventually remove them 
in the latter, as it makes A(k) less anisotropic. 

In this Letter, we investigate the pairing symmetry of 
KFe 2 As2 using thermal conductivity, a bulk directional 
probe of the superconducting gap [17] . All aspects of heat 
transport are found to be in agreement with theoretical 



expectation for a (i-wave gap [T5] [JJ5] , and inconsistent 
with accidental line nodes, whether vertical or horizontal. 
Moreover, the critical scattering rate r c for suppressing 
T c to zero is of order T c0 , as expected for d-wave, while 
it is 50 times T c0 in optimally-doped BaFe2As2 [2D] . 

Methods. - Single crystals of KFe2As2 were grown from 
self flux |21j . Two samples were measured: one for cur- 
rents along the a axis, one for currents along the c axis. 
Their superconducting temperature, defined by the point 
of zero resistance, is T c = 3.80±0.05 K and 3.65±0.05 K, 
respectively. Since the contacts were soldered with a su- 
perconducting alloy, a small magnetic field of 0.05 T was 
applied to make the contacts normal and thus ensure 
good thermalization. For more information on sample 
geometry, contact technique and measurement protocol, 
see ref. [6j. 

Resistivity.- To study the effect of impurity scatter- 
ing in KFe2As2, we performed measurements on a sin- 
gle crystal whose residual resistivity ratio (RRR) is 10 
times larger than that of the sample studied by Dong et 
al. [7] (Fig. la). To remove the uncertainty associated 
with geometric factors, we normalize the data of Dong 
et al. to our value at T = 300 K. A power-law fit below 
16 K yields a residual resistivity p a — 0.21 ± 0.02 pVt cm 
(2.24 ± 0.05 cm) for our (their) sample, so that 
p(300 K)/p = 1180 and 110, respectively. 

We attribute the lower p in our sample to a lower con- 
centration of impurities or defects. Note that except for 
the different po, the two resistivity curves p{T) are essen- 
tially identical (Fig. lb). Supporting evidence for a dif- 
ference in impurity /defect concentration is the difference 
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FIG. 1: (a) Electrical resistivity of the two samples of 
KFe2As2 studied here, with J || a (full red circles, left axis) 
and J || c (full blue squares, right axis). Our a-axis data is 
compared to that of Dong et al. [7] (open circles, left axis), 
normalized here to have the same value at T = 300 K (see 
text). The lines are a fit to p = po + aT a from which we 
extrapolate po at T — 0. (b) Same data for the two a-axis 
samples, up to 300 K. (c) Abrikosov-Gorkov formula for the 
decrease of T c with scattering rate T (line), used to obtain a 
value of F/r c for the three samples of KFe2As2, given their T c 
values and the factor 10 in po between the two a-axis samples 
(circles), assuming a disorder-free value of T c n = 3.95 K. 



in critical temperature: T c = 3.80±0.05 K (2.45±0.10 K) 
for our (their) sample. Assuming that the impurity scat- 
tering rate r a po, we can use the Abrikosov-Gorkov 
formula for the drop in T c vs T to extract a value of T/r c 
for the two samples, where T c is the critical scattering 
rate needed to suppress T c to zero (Fig. lc). We get 
r/r c = 0.05 (0.5) for our (their) sample. 

The c-axis resistivity p c {T) has the same temperature 
dependence as p a (T) below T ~ 40 K (Fig. la), with an 
intrinsic anisotropy Ap c /Ap a — 25 ± 1, where Ap = 
p{T) — po, with pco — 13 ± 1 pVL cm. We attribute the 
larger anisotropy at T — > 0, p c o/p a o = 60 ±10, to a larger 
r in our c-axis sample, consistent with the lower value of 
T c , from which we deduce r/r c = 0.1 (Fig. lc). 

Universal heat transport.- The thermal conductivity 
is shown in Fig. 2. The residual linear term kq/T is 
obtained from a fit to k/T = a + bT a below 0.3 K, 
where a = kq/T. The dependence of kq/T on magnetic 
field H is shown in Fig. 3. Extrapolation to H = 
yields n a0 /T = 3.6 ± 0.5 mW/K 2 cm and k c0 /T = 
0.18 ± 0.03 mW/K 2 cm. We compare to the data by 
Dong et al. [7j , normalized by the same factor as for elec- 
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FIG. 2: Thermal conductivity of KFe2As2, plotted as k/T 
vs T 2 , for J || a (« a , circles, left axis) and J || c (k c , squares, 
right axis), for a magnetic field H || c as indicated. Our a- 
axis data is compared to that of Dong et al. [7] (open circles, 
left axis), normalized by the same factor as in Fig. 1 (see 
text). Lines are a fit to n/T — a + bT a , used to extrapolate 
the residual linear term a = Kq/T at T = 0. For our a-axis 
sample (full red circles), a — 2.0, while for others a < 2. 



trical transport, giving Ko a /T = 3.32 ±0.03 mW/K 2 cm. 
At H — > 0, K-ao/T is the- same in the two samples (inset 
of Fig. 3), within error bars. 

This universal heat transport, whereby Ko/T is inde- 
pendent of the impurity scattering rate, is a classic signa- 
ture of line nodes imposed by symmetry |18[ I19j . Calcu- 
lations show the residual linear term to be independent 
of scattering rate and phase shift [18: , and free of Fermi- 
liquid and vertex corrections |19j . For a quasi-2D <i-wave 
superconductor [TB"l[T9"j : 
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where 7n is the residual linear term in the normal-state 
electronic specific heat, v-p is the Fermi velocity, and the 
superconducting gap A = A o cos(20) [2"2"] . 

ARPES measurements on KFe2As2 reveal a Fermi sur- 
face with three concentric hole-like cylinders centered 
on the r point of the Brillouin zone, labeled a, (3 
and 7, and a 4th cylinder near the X point [23l [24] . 
dHvA measurements detect all of these surfaces except 
the /3, and obtain Fermi velocities in reasonable agree- 
ment with ARPES dispersions, with an average value 
of vf — 4 x 10 6 cm/s The measured effective 

masses account for approximately 80% of the measured 
7n = 85 ± 10 mJ/K 2 mol [26l [27]. In d-wave symmetry, 
the gap in KFe2As2 will necessarily have nodes on all T- 
centered Fermi surfaces, and possibly on the X-centered 
surface as well [I4"| . The total kq/T may be estimated 
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FIG. 3: Field dependence of Ko/T obtained as in Fig. 2 (with 
corresponding symbols). The dashed line is a theoretical cal- 
culation for a d-wave superconductor with hT /Ao = 0.1 [38] . 
Inset: Zoom at low field. Lines are a power-law fit to extract 
the value of kq/T at H = 0. 



from Eq. 1 by using the average vp and the measured 
(total) 7n, which yields k qq /T = 3.3 ± 0.5 mW/K 2 cm, 
assuming A = 2.14 k B T c0 , with T c0 = 3.95 K. This 
is in excellent agreement with the experimental value of 
k /T = 3.6 ± 0.5 mW/K 2 cm. 

To compare with cuprates, the archetypal d-wave 
superconductors, we use Eq. 1 expressed directly in 
terms of va, the slope of the gap at the node, namely 
kqq/T ~ (k B /3hc)(v-p/v/^), with c the interlayer separa- 



tion [HI [19]. The ratio vf/va was measured by ARPES 
on Ba 2 Sr 2 CaCu 2 0g +( 5 [2S], giving Vf/va — 16 at opti- 
mal doping, so that k 00 /T ~ 0.16 mW/K 2 cm. This 
is in excellent agreement with the experimental value of 
k /T = 0.15±0.01 mW/K 2 cm measured in YBa 2 Cu 3 0. y 
at optimal doping [29]. 

In Fig. 4a, we plot k /T vs T for both KFe 2 As 2 
and YBa 2 Cu307, the superconductor in which univer- 
sal heat transport was first demonstrated [3D] . We see 
that Ko/T remains approximately constant up to at least 
hT ~ 0.5 fceTco in both cases. We conclude that both 
the magnitude of kq/T in KFe 2 As 2 and its insensitivity 
to impurity scattering are precisely those expected of a 
d-wave superconductor. By contrast, in an extended s- 
wave superconductor, there is no direct relation between 
kq/T and Ao, and a strong non-monotonic dependence 
on r is expected, since impurity scattering will inevitably 
make Ao less anisotropic [16]. This is confirmed by cal- 
culations applied to pnictides, which typically find that 
kq/T vs r first rises, and then plummets to zero when 
nodes are lifted by strong scattering [31] (see Fig. 4a). 

Critical scattering rate.- In a d- wave superconduc- 
tor, the critical scattering rate T c is such that hT c ~ 
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FIG. 4: Dependence of Ko/T (a) and T c (b) on impurity 
scattering rate F, normalized by T c o, the disorder-free super- 
conducting temperature, (a) Ko/T for KFe2As2 (red circles; 
see text) and the cuprate YBa2Cu307 (blue squares; from 
ref. !3Up. normalized by the theoretically expected value for a 
d-wave superconductor, koo/T — 3.3 and 0.16 mW/K 2 cm, re- 
spectively (see text). The typical dependence expected of an 
s-wave state with accidental nodes is also shown, from a cal- 
culation applied to pnictides (black line; from ref. I31[). (b) T c 
for KFe2As2 (red circles; from Fig. lc) and for the pnictides 
BaFe2As2 and SrFe2As2 at optimal doping (from ref. [20]) . 



^bTco [32 . We can estimate T c for KFe 2 As 2 from the 
critical value of po, evaluated as twice that for which 
r/r c = 0.5 in Fig. lc, namely pg 1 ~ 4.5 p£L cm. Us- 
ing Lo/p"^ = 7n«f t c/ 3 > where L = (tt 2 / '3) (fc B / e) 2 , we 
get HT C = h/2r c ~ 1.3 ± 0.2 ksTcQ, in excellent agree- 
ment with expectation for a d-wave state. By contrast, 
hV c /kBTco — 45 in BaFe 2 As 2 and SrFe 2 As 2 at optimal 
Co, Pt or Ru doping [H] (see Fig. 4b). This factor 30 
difference in the sensitivity of T c to impurity scattering is 
proof that the pairing symmetry of KFe 2 As 2 is different 
from the s-wave symmetry of Co-doped BaFe 2 As 2 [B]. 

Direction dependence. In the case of a <i-wave gap on 
a single quasi-2D cylindrical Fermi surface (at the zone 
center), the gap would necessarily have 4 line nodes 
that run vertically along the c axis. In such a nodal 
structure, zero-energy nodal quasiparticles will conduct 
heat not only in the plane, but also along the c axis, 
by an amount proportional to the c-axis dispersion of 
the Fermi surface. In the simplest case, c-axis conduc- 
tion will be smaller than a-axis conduction by a fac- 
tor equal to the mass tensor anisotropy (vp in Eq. 1). 
In other words, {n a0 /T) / '(k c0 /T) ~ {k^/T)/ (k c n/T) = 
(caN)/(ccN), the anisotropy in the normal-state thermal 
and electrical conductivities, respectively. This is con- 
firmed by calculations for a quasi-2D d-wave supercon- 
ductor [34] , whose vertical line nodes yield an anisotropy 
ratio in the superconducting state very similar to that 
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of the normal state. This is what we see in KFe2As2 
(inset of Fig. 3): {n aQ /T)/{n cQ /T) = 20 ± 4, very close 
to the intrinsic normal-state anisotropy (c a N)/(o'cN) = 
(Ap c ) / (Ap a ) = 25 ± 1. This shows that the nodes in 
KFe2As2 are vertical lines running along the c axis, rul- 
ing out proposals [TT] of horizontal line nodes lying in a 
plane normal to the c axis. 

Moreover, the fact that the Fermi surface of KFe2As2 
contains several sheets with very different c-axis disper- 
sions 25 , 35 provides compelling evidence in favor of d- 
wave symmetry. In an extended s-wave scenario, the gap 
would typically develop vertical line nodes on some but 
not all zone-centered sheets of the Fermi surface [15] , and 
so the anisotropy in k would typically be very different 
in the superconducting and normal states, unlike what 
is measured. By contrast, in d-wave symmetry all zone- 
centered sheets must necessarily have nodes, thereby en- 
suring automatically that transport anisotropy remains 
similar in the superconducting and normal states. 

Temperature dependence.- So far, we have discussed 
the limit T — > and H — > 0, where nodal quasiparticles 
are excited only by the pair-breaking effect of impuri- 
ties. Raising the temperature will further excite nodal 
quasiparticles. Calculations for a d-wave superconductor 
show that the electronic thermal conductivity grows as 

T 2 pil[22]: 



where a is a dimensionless number and is the impurity 
bandwidth, which grows with the scattering rate T [18]. 
A T 2 slope in k/T was resolved in YBa2Cu307 [29]. 

Our KFe2As2 sample shows a clear T 2 dependence be- 
low T ~ 0.3 K, with K a /T = (n a0 /T){l + 23 T 2 ) (Fig. 2). 
Comparison with the data by Dong et al. [7] reveals that 
this T 2 term must be due to quasiparticles. Indeed, be- 
cause phonon conduction at sub-Kelvin temperatures is 
limited by sample boundaries and not impurities |33j , the 
fact that the slope of k/T in their sample (of similar di- 
mensions) is at least 10 times smaller (Fig. 2), implies 
that the larger slope in our data must be electronic. 

In the limit of unitary scattering, 7 2 oc T, so that 
a 10-times larger T would yield a 10-times smaller T 2 
slope [18] , consistent with observation. The tempera- 
ture below which the T 2 dependence of n e /T sets in, 
T ~ 0.1 T c , is a measure of 7. It is in excellent agree- 
ment with the temperature below which the penetration 
depth A a (T) of KFe2As2 (in a sample with similar RRR) 
deviates from its linear T dependence [5] , as expected of a 
d-w&ve superconductor [35] . Note that the T dependence 
of k/T for an extended s-wave gap is not T 2 [3T] . 

Magnetic field dependence.- Increasing the magnetic 
field is another way to excite quasiparticles. If the gap 
has nodes, the field will cause an immediate rise in 
kq/T [TT] |371 [3B], as observed in all three samples of 
KFe2As2 (inset of Fig. 3). Calculations for a d-wave 



superconductor in the clean limit (hT <C k^T c ) yield a 
non-monotonic increase of kq /T vs H [38] in remarkable 
agreement with data on the clean sample (Fig. 3). 

A rapid initial rise in kq/T vs H has been observed 
in the cuprate superconductors YBa2Cu307 [31] and 
Tl 2 Ba 2 Cu0 6+5 [D]. In the dirty limit, KFe 2 As 2 [7J and 
Tl 2 Ba2Cu06+5 [40] show nearly identical curves of k /T 
vs H/H c2 (see ref. [7J. Measurements on cuprates in the 
clean limit, such as optimally-doped YBa2Cu30 !/ , have 
so far been limited to H <C H C 2- 

In summary, all aspects of the thermal conductivity of 
KFe 2 As 2 , including its dependence on impurity scatter- 
ing, current direction, temperature and magnetic field, 
are in detailed and quantitative agreement with theoret- 
ical calculations for a d-wave superconductor. The scat- 
tering rate needed to suppress T c to zero is exactly as 
expected of d-wave symmetry, and vastly smaller than 
that found in other pnictide superconductors where the 
pairing symmetry is believed to be s-wave. This is com- 
pelling evidence that the pairing symmetry in this iron- 
arsenide superconductor is d-wave, in agreement with 
renormalization-group calculations |14j . Replacing K in 
KFe2As2 by Ba leads to a superconducting state with a 
10 times higher T c , but with a full gap without nodes [3], 
necessarily of a different symmetry. Understanding the 
relation between this factor 10 and the pairing symmetry 
provides insight into what boosts T c in these systems. 
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